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Abstract. We review recent results on the production of prompt charmonium in association with a 
hadron jet or a prompt photon in two-photon collisions at next-to-leading order in the factorization 
framework of nonrelativistic quantum chromodynamics. 
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^ : 1. INTRODUCTION 

o i 

The factorization formalism of nonrelativistic quantum chromodynamics (NRQCD) [1] 
provides a rigorous theoretical framework for the description of heavy-quarkonium 
production and decay that is renormalizable and predictive. Theoretical predictions 
are decomposed into sums over products of short-distance coefficients, which can be 
calculated perturbatively as expansions in the strong-coupling constant a s , and long- 
distance matrix elements (MEs), which are subject to relative- velocity (v) scaling rules 
Ph. and must be extracted from experiment, and are so organized as double expansions in 

a s and v. This formalism takes into account the complete structure of the QQ Fock 



X 



space, which is spanned by the states n = 2S + l Ly> with definite spin S, orbital angular 
momentum L, total angular momentum /, and color multiplicity a = 1,8, and so predicts 
the existence of color-octet (CO) processes in nature. 

The greatest triumph of the NRQCD factorization formalism was its ability to cor- 
rectly describe the cross section of inclusive charmonium hadroproduction at the Teva- 
tron, which exceeds the color-singlet-model prediction by more than one order of mag- 
nitude. In order to convincingly establish the phenomenological significance of the CO 
processes, it is indispensable to identify them in other kinds of high-energy experiments 
as well. The verification of the NRQCD factorization hypothesis is presently hampered 
both from the theoretical and experimental sides. On the one hand, the theoretical pre- 
dictions to be compared with existing experimental data are, apart from very few excep- 
tions, of lowest order (LO) and thus suffer from considerable uncertainties. The mea- 
surement of charmonium polarization at the Tevatron currently presents a challenge for 
NRQCD factorization, but any conclusions are premature in the absence of a full-fledged 
next-to-leading-order (NLO) analysis. It is, therefore, mandatory to calculate the NLO 
corrections to the hard- scattering cross sections and to include the composite operators 
that are suppressed by higher powers in v. Apart from the usual reduction of the renor- 



malization and factorization scale dependences, sizeable effects, e.g. due to the opening 
of new partonic production channels, are expected at NLO. On the other hand, the ex- 
perimental errors are still rather sizeable. The latter are being significantly reduced by 
HERA II and run II at the Tevatron, and will be dramatically more so by the LHC and 
hopefully a future e + e~ linear collider such as the TeV-Energy Superconducting Linear 
Accelerator (TESLA), which is presently being designed and planned at DESY. 

Recently, 2^2 processes of heavy-quarkonium production were for the first time 
studied at NLO in the NRQCD factorization formalism [2, 3]. Specifically, the produc- 
tion of prompt charmonium, which is produced either directly or through the decay of 
heavier charmonia, with finite transverse momentum (pj) in association with a hadron 
jet [2] or a prompt photon [3] via direct photoproduction in two-photon collisions was 
considered. In this presentation, we review the most important conceptional issues and 
phenomenological results of Refs. [2, 3], in Sections 2 and 3, respectively. 

2. CONCEPTIONAL ISSUES 

We focus attention on the process yy — > where J / y is promptly produced at 

finite value of p T and X is a purely hadronic remainder. Since the incoming photons can 
interact either directly with the quarks participating in the hard- scattering process (direct 
photoproduction) or via their quark and gluon content (resolved photoproduction), this 
process receives contributions from the direct, single-resolved, and double-resolved 
channels, which are formally of the same order in the perturbative expansion. At LO, 
the bulk of the cross section is due to single-resolved photoproduction, and the NRQCD 
prediction [4] based on the MEs determined from fits [4] to Tevatron data nicely agrees 
with a recent measurement by the DELPHI Collaboration at LEP2 [6]. 

Here, we consider direct photoproduction at NLO [2]. At LO, there is only one 
partonic subprocess, namely 

y+y^ccfsf^+g. (1) 

At NLO, virtual corrections to process (1) and (0 H [ 3 sf^]), where H = J/y^Xcj, V > and 
real corrections to 

y + y - cc[n] +g + g, n = ¥/\ tf 8) , 3 SS 8) , 3 /f , (2) 

Y+Y^cc[n]+u + u, n = 3 sf\ (3) 

Y + Y^cc[n)+q + q, n= 1 S { * ) , 3 S?\ 3 PP, (4) 

where u and u denote the Faddeev-Popov ghosts of the gluon, contribute. 

The virtual corrections to process (1) receive contributions from self-energy, trian- 
gle, box, and pentagon diagrams. The self-energy and triangle diagrams are in general 
ultraviolet (UV) divergent; the triangle, box, and pentagon diagrams are in general in- 
frared (IR) divergent; and the pentagon diagrams without three-gluon vertex also con- 
tain Coulomb divergences. As for the light-quark loops, the triangle diagrams vanish by 
Furry's theorem, while the box diagrams form a finite subset. The virtual corrections to 



(<ff H [ 3 S\ ]) also produce UV, IR, and Coulomb divergences. The UV and IR divergences 
are extracted using dimensional regularization in d = 4 — 2e space-time dimensions, 
leading to poles in £uv and £ir, respectively, while the Coulomb singularities are regu- 
larized by a small value of v. The UV divergences are removed by the renormalization of 

{& H [ 3 S^]), a s , the charm-quark mass and field, and the gluon field, which is performed 
in the modified minimal- subtraction (MS) scheme for the former two quantities, render- 
ing them dependent on the renormalization scales A and /i, respectively, and in the on- 
mass-shell scheme for the residual three quantities. The IR divergences cancel among 

the virtual and real corrections, the wave-function renormalizations, and (<f? H [ 3 S^]). 

The Coulomb divergences cancel between the virtual corrections and (€? H [ 3 S^]}. 

The real corrections are plagued by IR divergences, which come as collinear diver- 
gences from the initial state and collinear and/or soft ones from the final state. They are 
identified by appropriately slicing the three-particle phase space using small parameters 
8i and 8f, respectively. The collinear and/or soft regions of phase space are integrated 
over analytically in d dimensions, while the hard region is integrated over numerically 
in four dimensions. The sum of these contributions is, to very good approximation, in- 
dependent of 8i and 8f. The initial-state collinear divergences are factorized at some 
factorization scale M and absorbed into the parton density functions (PDFs) of the q 
and q quarks inside the resolved photon. The M dependence thus introduced is approxi- 
mately compensated by the LO single-resolved contribution. 

Combining the contributions arising from the virtual corrections (vi), the parameter 
and wave-function renormalization (ct), the operator redefinition (op), the initial-state 
(is) and final-state (fs) collinear configurations, the soft-gluon radiation (so), and the 
hard-parton emission (ha) as 

do O , A , M) = do ( fi , A ) [ 1 + Svi (ju ; £uv, £ir, v) + 5 ct (ju ; euv , £ir) + 8 op ( ju , A ; em., v) 
+ 5f S (ju; £ir, 8 f )] + do^H, A, Af ; $) + do so (n, A; £ir, 8 f ) + da h& (fi,X;8i, 8 f ) , (5) 

the regulators £uv, £ir, v, <5,, and Sf drop out and the /i and A dependences formally 
cancel up to terms beyond NLO, while the M dependence is unscreened at NLO. 

3. PHENOMENOLOGICAL RESULTS 

We consider two-photon collisions at TESLA operating at a center-of-mass energy of 
500 GeV, where the photons arise from electromagnetic initial-state bremsstrahlung, 
with antitagging angle max = 25 mrad, and beamstrahlung, with effective beam- 
strahlung parameter T = 0.053. The J fy, % c j, and y/ 1 MEs are adopted from Ref. [5] 
and the photon PDFs from Ref. [7]. 

In Fig. 1, we study d 2 o /dpjdy (a) for rapidity y = as a function of pj and 
(b) for pr = 5 GeV as a function of y, comparing the LO (dashed lines) and NLO 
(solid lines) results of direct photoproduction with the LO result of single-resolved 
photoproduction (dotted lines). From Fig. 1(a), we observe that, with increasing value 
of pj, the NLO result of direct photoproduction falls of considerably more slowly 
than the LO one. In fact, the QCD correction (K) factor, defined as the NLO to LO 



ratio, rapidly increases with pj, exceeding 10 for /T^IO GeV. This feature may be 
understood by observing that so-called fragmentation-prone [8] partonic subprocesses 
start to contribute to direct photoproduction at NLO, while they are absent at LO. Such 
subprocesses contain a gluon with small virtuality, q 2 = Am 2 , that splits into a cc pair in 

the Fock state n = 3 S^ and thus generally generate dominant contributions at pj ^> 2m c 
due to the presence of a large gluon propagator. In single-resolved photoproduction, a 
fragmentation-prone partonic subprocess already contributes at LO. This explains why 
the solid and dotted curves in Fig. 1(a) run parallel in the upper pj range. At low values 
of pt, the fragmentation-prone partonic subprocesses do not matter, and the relative 
suppression of direct photoproduction is due to the fact that, at LO, this is a pure CO 
process. 
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FIGURE 1. LO single-resolved (dotted lines), LO direct (dashed lines), and NLO direct (solid lines) 
contributions to d 2 jdpjdy (a) for y — as a function of pr and (b) for pj — 5 GeV as a function of y. 

In the case of 77 — > //i/z + Xy, where Xy contains a prompt photon, the K factor was 
found to decrease fast with increasing value of p T , falling below 01 for p T ^ 14 GeV [3]. 
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